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ABSTRACT 

In solar active regions (ARs), Doppler shifts measured along polarity inversion lines (PILs) of the line- 
of-sight (LOS) magnetic field determine one component of the velocity perpendicular to the magnetic 
field. Along PILs, these velocities can be used to: (i) improve estimates of electric fields, which can 
be used to derive the Poynting flux of magnetic energy across the photosphere; and (ii) constrain the 
physical processes underlying flux cancellation, the mutual apparent loss of magnetic flux in closely 
spaced, opposite-polarity magnetogram features. At least two factors make the zero point of measured 
Doppler velocities uncertain: instrumental variations (from, e.g., thermal effects); and the convective 
blueshift, a known correlation between intensity and upflows. We present a method to absolutely 
calibrate LOS velocities using three successive vector magnetograms and one Dopplergram coincident 
with the central magnetogram. It assumes ideal electric fields govern magnetic field evolution along 
PILs, and optimizes consistency between changes in LOS flux near PILs and the transport of transverse 
magnetic flux by LOS velocities. Our method can account for residual convective blueshift at the 
limb, while an alternative method, fitting the center-to-limb variation of Doppler velocities, cannot. 
Applying the method to vector magnetograms of AR 11158 from the Helioseismic and Magnetic Imager 
(HMI) on the Solar Dynamics Observatory, observed by the Helioseismic and Magnetic Imager (HMI) 
aboard the Solar Dynamics Observatory, we find clear evidence of offsets in the Doppler zero point, 
in the range of 50 - 550 m s -1 . In principle, our method can also be used to detect non-ideal electric 
fields. We also discuss shortcomings in our initial implementation, and suggest ways to address these. 
In addition, we note that a simpler calibration can be determined from an LOS magnetogram and 
Dopplergram pair from the median Doppler velocity among all near-disk-center PIL pixels. 



1. INTRODUCTION 

Solar variability is intimately related to magnetic flux 
at the solar photosphere: solar flares, coronal mass ejec- 
tions (CMEs), and enhanced radiation from solar activ- 
ity, ranging from radio to X-ray wavelengths, all occur in 
the outer solar atmosphere above magnetized regions of 
the photosphere. Fundamentally, these manifestations of 
solar activity are driven by the release of energy stored 
in or transmitted by magnetic fields. 

1.1. Doppler Shifts & the Poynting Flux 

Magnetic energy passes from the solar interior into the 
Sun's outer atmosphere as an outward-directed Poynting 
flux S of magnetic energy across the photosphere, 

S n = n- S = cn- (E x B)/4tt , (1) 

where E and B are the photosphcric electric and mag- 
netic fields, and n is the outward-directed unit vector 
normal to the photosphere. We will refer to vectors per- 
pendicular to the local normal as horizontal. We use the 
photosphere as the boundary between the solar interior 
and outer atmosphere primarily for convenience, because 
the vector magnetic field is only routinely measured at 
the photosphere. A different atmospheric layer could be 



used if observations of the vecto r magnetic field were rou- 
tinely produced there (see, e.g., Mctcalf et al.lll995|) . 

While B at the photosphere can be measured by vec- 
tor magnetographs, E must be inferred by other means. 
A popular approach has been to use observed magnetic 
evolution in sequences of magnetograms to estimate the 
electric field by inverting the finite-difference approxima- 
tion to Faraday's law, 

AB 

— = -c(V x E) , (2) 

where At is the difference between the times tf and U of 
final and initial magnetograms, respectively, with E rep- 
resenting an average electric field over At. Note that for 
any E which satisfies this expression, (E — Vip), where ip 
is an arbitrary scalar potential function, will also satisfy 
it. Hence, Faraday's law does not fully constrain E. 

Assuming the magnetic field is frozen to the plasma at 
the photosphe re (for argum ents supporting this assump- 
tion, see, e.g., lParkerll!984f) . the ideal Ohm's law relates 
photospheric velocities, v, to the electric field, 

cE = -v x B . (3) 

This implies estimates of v can be used to deter- 
mine the flux of magnetic energy (and magnetic helic- 
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ity) across the pho tosphere (e.g., iDemoulin and Bergerl 
l2003t iSchuckl [2006) . Several techniques have been de- 
velop ed to estima t e photospheric flows fr om AB n /At, 
e.g.. [Chad (120011): iKusano et al.1 (120021): IWelsch et alJ 
(120041): I Schuckl (120061 ): iFisher and Welschl (j2008fl and 
ISchuckl (l2008fl . These techniques are, however, imperfec t 
(IRieutord etaH[200lt IWelsch et all [20071: ISchuckl 1200^ . 
s o efforts to improve them are ongoing. 

IFisher et al.l (|2010[) presented a method to determine E 
from a sequence of vector magnetograms using a poloidal- 
toroidal decomposition (PTD) of the magnetic field, with 
Faraday's law. Notably, the PTD method uses evolution 
of both the normal magnetic field, AB n / At, and the 
normal electric current, AJ n / At, to estimate E. While 
the PTD approach does not rely upon the ideal Ohm's 
law to estimate E, any Ohm's law with known resistive 
terms, including the ideal case, can be imposed post facto 
to constrain the soluti on for E. 

For the ideal case, IFisher et al.l (|2012[ ) recently pre- 
sented a method to use Doppler measurements of the 
line-of-sight (LOS) component of the velocity, vlos, to 
better constrain E. They tested this approach using syn- 
thetic magnetograms and Doppler data from an MHD 
simulation. They noted that flows along B (commonly 
referred to as parallel or siphon flows) do not contribute 
to the ideal electric field in equation ([3]), but can con- 
tribute to vlos m regions where the LOS component of 
the magnetic field is nonzero (Blos 7^ 0). Hence, they 
only incorporated Doppler data from areas: (i) near po- 
larity inversion lines (PILs), loci where the component 
of the magnetic field along their (assumed) LOS changes 
sign and the LOS field vanishes; and (ii) where the field 
transverse to the LOS is large. Their tests demonstrated 
that including Doppler data near PILs substantially im- 
proves estimation of both E and the normal Poynting 
flux S n . For ideal evolution, this makes sense because, 
in principle, the Doppler electric field ~E D from the LOS 
velocity i>los an d transverse magnetic field B trs along a 
PIL, 
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-(fLOS X B tra ) | PIL > 



(4) 



is not uncertain by the gradient of a scalar potential, as 
are estimates of E from equation © alone. 

In real magnetograms, procedures have been devel- 
oped to automatically i dentify PILs (also s ometimes 
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Falco ner et alj 20031) of both 
Schriiverll20'07tlWelsch and Lil 

and normal components 



called neutral lines; e.g., 
the LO S components (e.g., 
2001 IWelsch et al.l 120091) 
Falconer et al.ll2003l) of magnetogram fields 

One factor hampering studies of Poynting fluxes 
has been the relative dearth of sequences of vector 
magnetograms. The launch of the SpectroPolarime- 
ter (S P) instrument with the Solar Optical Telescope 
(SOT: iTsuneta, et al.1 [2001 ) aboard the Hinode satellite 
(jKosugi et al.l 120071) provided some seeing-free magne- 
togram sequences for investigations of PIL dynamics, 
such as obse rvations suggest i ve of a n emerging flux rope 
reported by lOkamoto et al.l (j2008f l. However, the ca- 
dence of the SP instrument is relatively slow compared to 
timescales of photospheric evolution on arcsecond scales, 
and the field of view (FOV) of SOT is limited. Vector 
magnetograms of active region fields at higher cadence 
and over larger FOVs should be routinely produced us- 
ing the Helioseismic and Magnetic Imager (HMI) instru- 



ment (Schcrrc r et al.l [20121 ) aboard the Solar Dynamics 
Observatory (SDO), enabling routine estimates of Poynt- 
ing flux. 

1.2. Emergence & Cancellation Along PILs 

Given the association between solar activity and mag- 
netic fields, a quantitative description of processes re- 
sponsible for the introduction and removal of magnetic 
flux into the solar atmosphere is essential to understand 
solar activity. 

New magnetic flux typically appears at the solar pho- 
tosphere via emergence of il-shaped loops from the con- 
vection zone (Figured! top row, left panel). About 3000 
active regions typically emerge over one 11-year cycle of 
solar activity, each with ~ 10 22 Mx (10 14 Wb) in un- 
signed flux. In the quiet Sun, estimates vary for the time 
required for emergence to replace the small-scale, tur- 
bulent magnetic fields there, which totals ~ 10 23 Mx in 
unsigned flux over the photospheric surface; but they are 
generally on the order of a day or less (e.g., Hagenaar et 
al. 2003). We note that new flux can also appear in 
photospheric magnetograms via the submergence of the 
bases of U-shap e d loop s from above, as seen in simula- 
tions bv lAbbettl ([2001 (Figure Q] top row, right panel), 
though this process is not often discussed. Hence, the 
term appearance encompasses both emergence and sub- 
mergence. This definitio n of the term diffe rs substan- 
tially from that used by lLamb et al.l ([2008), in which 
appearance refers to initial identification of a feature by 
a tracking algorithm, where the definition of "feature" is 
algorithm-dependent. 

Much like the appearance of new flux, the removal of 
magnetic flux from the photosphere also plays a central 
role in solar activity. Clearly, for the large-scale solar 
dynamo to operate cyclically, flux that is introduced to 
the photosphere must eventually be removed. Similarly, 
flux emergence in the small-scale dynamo must be stat is- 
tically balanced by flux removal ([Schriiver et al.|[l997D . 

How is flux removed from the photo sphere? Th e short 
answer is "flux cancellation," which iLivi et all ([19851 ) 
defined in observational terms as "the mutual appar- 
ent loss of magnetic flux in closely spaced features of 
opposite polarity" in magnetogram sequences. Phys- 
ically, cancellation could correspond to (i) the emer- 
gence of U-shaped magnetic loops (e.g., iLites et ail 
119951 : Ivan Driel-Gesztelvi et all 120001). Qi) the sub- 
mergence of f^-shaped loops (e.g.. [Rabin et all Il984t 
IHarvev et all [19991 IChae et al.l 1 120041 : llida et al.1 12010F . 
or (iii) reconnection in the magnetogram layer (e.g. , 
I Yurchvshvn and Wand 120011: iKubo and Shimizul 120071: 
see also IWelsch! 120061) . These possibilities are sketched 
i n the bottom r ow of Figure Q] 

ISpruit et al.l (jl987| ) proposed that opposite polarities 
within active regions could be connected by U-loop sub- 
surface extensions, and that flux in these extensions 
might emerge as weak, "sea-serpent" (undulating) fields 
between strong field regions, eventually ca nceling with 
active region flux via U-loop emergence. Eo^HqqiI) 
proposed a process that also removes most active region 
flux from the photosphere by the emergence of U-loops. 
He suggests that the flux rope that forms an active re- 
gion contains many turns, such that a magnetogram can 
intersect the flux rope many times after the flux rope 
has partially emerged, and as the flux rope continues 



Figure 1. Top row: sketch of field lines (thin, solid) in two possible processes that cause more flux to appear in magnetograms, U-loop 
emergence (left) and U-loop submergence (right), viewed in a plane perpendicular to the polarity inversion line (PIL; vertical, dashed line 
in all panels). The horizontal solid line represents the magnetogram layer. Near disk center, if the evolution is ideal, then the plasma 
velocity (fat, white arrows) should produce Doppler shifts along the PIL. Bottom row: Sketches of three possible processes in magnetic flux 
cancellation: emergence of U-loops (left), submergence of fl-loops (middle), or reconnective cancellation (right). As with the top row, if the 
evolution is ideal, then the plasma velocity (fat, white arrows) should produce Doppler shifts along the PIL for emergence or submergence 
flux (left and middle panels, respectively) near disk center. In principle, processes that increase or decrease total unsigned flux only occur 
along PILs. 



to emerge, much of the flux cancels by U-loop emer- 
gence, with successive magnetograms showing less and 
less flux rem aining until little is left. In contra st to 
these models. Ivan Ballegooiien and Mackavl (|2007[ ) and 
Ivan Ballegooiie n (2008) constructed numerical models of 
active region flux systems remaining anchored near the 
base of the convection zone, with most canceling flux re- 
tracting to "repair" the toroidal magnetic field deep in 

t he solar interi or. 

iSpruit et al.1 (|1987l ) noted that sea-serpent cancel- 
lation might occur on small scales. On unobserv- 

ably small scales, this would produce appare nt in 

situ flux disappearance (IWallenhorst and Howard) |1982l ; 
IWallenhorst and Topkal Il982t IGaizauskas et al.l I1983T ). 
Further, unlike t he "self-canc e llatio n" of active region 
flux modeled bv ISpruit et all (119871) , lLowl 1)200 ID . and 
van Ballegooiien and Mackay (2007), active region fields 
might also cancel with the ubiquitous, small-scale fields 
of the quiet sun (e.g., Lin & Rimmel e 1999, Harvey et 
al. 2 007), including unresolved fields ([Sanchez Almeidal 
2009). This process would formally remove equal 
amounts of active-region and quiet-Sun flux from the 
photosphere, but if the latter were below a given mag- 
netograph's noise threshold, it would be undetectable, 
and the active-region flux would have seemed to dis- 
appear. While case studies of cancellation have been 
made, it is unclear how much active region flux cancels 
on observable scales over a solar cycle. Hence, it is pos- 
sible that cancellation with quiet-Sun flux is the domi- 
nant xnethod_of active region flux removal. We remark 
that lKubo et a l. (2008) noted that moving magnetic fea- 
tures in moat flows around sunspots can explain rates 
of sunspot flux loss, but much flux in active regions lies 
outside sunspots, and its disappearance must still be un- 



derstood. 

Because magnetic fields are divergence- free, all field 
lines (tangent lines of the vector field B) form closed 
loops. (These might be infinitely long, or ergodic, but 
in any case field lines do not end.) Magnetic flux can 
therefore only emerge or cancel where magnetic fields 
are tangent to the photosphere. These loci correspond 
to normal-field PILs, where the normal magnetic field 
vanishes and regions of positive and negative flux, where 
emerging (or submerging) field lines thread the photo- 
sphere, are nearby. In our terminology, the appearance 
(or cancellation) of magnetic flux increases (resp., de- 
creases) the total unsigned magnetic flux at the photo- 
sphere. 

If the photosphcric electric field during cancellation is 
ideal or nearly so, then measurements of time-averaged 
Doppler shifts along PILs should be able to distinguish 
between cancellation via either U-loop emergence or Q- 
loop submergence. Lack of a clear Doppler signal while 
LOS flux cancels would be consistent with reconnective 
cancellation. We note that the magnetic field along PILs 
away from disk center in LOS magnetograms can have 
a component that is normal to the photosphere, im- 
plying LOS PILs away from disk center do not exactly 
correspond to sites of flux appearance (or cancellation) . 
Hence, only Doppler shifts along LOS PILs near disk 
center can effectively constrain the physical processes at 
work in cancellation. 

Several case studies of Dopp ler shifts at cancella- 
tion s ites hav e been undertaken. [Yurchyshyn and Wand 
(|2001l ) and iBellot Rubio and Beck! (|2005l ) reported 
Doppler shifts consistent with upflows at the cancellation 
sites they studied, which they interpreted as out flows 
from reconnective cancellation. iChae et al.l ()2004l ) and 
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llida et all (|2010l ) fou nd evidence for flu x subm ergence 
during cancellation. iKubo and Shimizul (|2007l ) studied 
cancellations along several PILs in ASP magnetograms, 
and reported mixed Doppler signals, which they inter- 
preted in terms of reconnective cancellation at multiple 
heights. To determine th e rest wavelengths used to com- 
pute Doppler velocities, iChae et ail (120041) used quiet - 
Sun values of line cen ter, while IKubo and Shimizul ()2007f ) 
and llida et al.l (|2010D estimated their rest wavelengths by 
averaging line centers over their FOVs. As described be- 
low, however, these approaches to determining rest wave- 
lengths probably yield biased estimates. 

1.3. Biases in Doppler Shifts Along PILs 

One factor complicating use of Doppler data both to 
estimate photospheric electric fields and to constrain dy- 
namics in flux cancellation is inaccuracies in determina- 
tion of the rest wavelengths of photospheric lines in ac- 
tive regions. Uncertainties in rest wavelength can arise 
from both instrumental effects (e.g., from thermal vari- 
ations in components) and biases in analysis techniques. 
Periodicities in magnetic fields estimated by HMI (Liu, 
Hoeksema, et al. 2012b) on orbital time scales (12 and 24 
hr) suggest instrumental effects are present in estimated 
magnetic fields. It is plausible that similar effects should 
be present in Doppler signals, a point we revisit below. 

Because the position of line center is typically 
computed from sampling predominantly quiet- 
Sun regions, where line profiles are systcmati- 
cally shifted blue-w a rd by the co n vectiv e blueshift 
(iDravins et al.l ll98lHCavallini et al.l 119861: IHathawavl 
119921: lAsplund and Colletl 120031; iSchuckl 120101 ) . Doppler 
shifts in active regions typically exhibit "pseudo- 
rcdshifts." The blue-ward bias of quiet-Sun line profiles 
occurs because rising plasma is both (i) brighter than 
sinking plasma (since rising plasma is hotter) and (ii) 
occupies a greater fraction of an instrument's FOV than 
sinking plasma (since upwelling convective cells are 
larger than downflow lanes). Consequently, a determi- 
nation of line center position based upon the statistical 
properties (e.g., means or medians) of Doppler images 
(Dopplergrams) is biased by the upward motion of quiet- 
Sun p lasma. Observations and mo deling , by, e.g.. iGravl 
(2009) and lAsplund and Colletl (|2003l i. respectively, 
suggest that the magnitude of this bias can range from 
a few hundred m sec" 1 to nearly 1 km sec -1 for various 
photospheric lines. Beca use active region ma gnetic fields 
inhibit convection (e.g., iWelsch et al.l 120121 ) . measured 
line centers in active regions are red-shifted relative to 
any rest wavelength derived from quiet-Sun Doppler 
measurements. (Helioseismology requires accurate mea- 
surement of changes in Doppler shifts, so is insensitive 
to errors in determination of the rest wavelength.) 
P. Scherrer (private communication 2009) has stated 
that one must account for the convective blueshift to 
accurately determine Doppler shifts in active regions. 

By fitting measured Doppler shifts over the disk with 
profiles that account for differential rotation, meridional 
flow, and center-to-limb variations, any overall constant 
Doppler s hift (sometimes referred to as the convective 
limb s hift: lHathawavlll992D can be est imated (jSnodgrass! 
fl98i IHathawavl 119921: iSchuckl 120101) . This approach, 
however, has a major physical uncertainty: the physics 
of the center-to-limb variation in average Doppler shift 



in the particular spectral line used by HMI (or any 
other spectral line) involves detailed interactions between 
height of formation, the height of convective turnover, 
the variation with viewing angle of the average convective 
flow speed , and the variation w ith viewing angle of opti- 
cal depth (jCarlsson et al.l 120041) . We note that diverging 
flows tangent to the photosphere in granules can, de- 
pending upon optical depth in granules at the formation 
height of the line, produce a convective blueshift toward 
the limb, because diverging flows on the near sides of 
granules approach the viewer, while receding flows on the 
far sides of the granules are at least partially obscured 
by the optical depth of the granules. This is sketched 
in Figure [2] Hence, fitting the observed center-to-limb 
variation in line-center positions does not imply that all 
bias from convective motions has been determined. In 
princi ple, simulations of convection near the photo sphere 
(e.g.. lAsplund and Colletl l200l iFleck et all 1201 If) could 
be used to characterize the expected center-to-limb vari- 
ation in line centers, which could then be used to remove 
the modeled convective blueshift bias (but not any in- 
strumental biases). We are unaware of published com- 
parisons of observed and modeled center-to-limb varia- 
tions in the HMI spectral line. 




Figure 2. Idealized convective motions near the photosphere from 
center to limb. Thick black arrows represent velocity vectors. Thin 
black solid lines represent stream lines. Thick dashed line repre- 
sents optical d epth unity, and shows granules' three-di mensional 
structure (e.g., Bcrgcr et al. 2007; Carlsson et al. 2004). For an 
observer looking from the top of the figure down, Doppler veloc- 
ities along lines of sight (thin vertical dashed lines ) both at disk 
center and toward the limb intersect more areas with velocities 
toward the observer than away. This is one reason that determina- 
tion of the rest wavelength by averaging (or taking the median ) of 
line center positions over the FOV will result in an estimate biased 
blue-ward of the true rest wavelength. This analysis ignores the 
additional factor that upBowing plasma is brighter than downBow- 
ing plasma, so the former contributes more to intensity-averaged 
line proBles. 

The goal of this paper is to describe a "magnetic cal- 
ibration" technique to estimate any spatially uniform 
bias in measured Doppler velocities. In basic terms, this 
method estimates such a bias by requiring statistical con- 
sistency between two independent measures of changes 
in flux near LOS-field PILs, the loci where flux emerges 
and submerges: (1) A$los, half the change of total un- 
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signed LOS flux near PILs over a time interval At; and 
(2) A& v bl, the amount of transverse flux transported 
upward or downward across the magnetogram layer over 
At, computed by summing the Doppler electric field — 
the product of the Doppler velocity and transverse mag- 
netic field strength — along PILs. This consistency relies 
upon both Faraday's Law and the assumption of ideal 
evolution. 

The rest of the paper is organized as follows. In 
the next section, we describe our calibration method in 
greater detail. We then apply the method to the ini- 
tial vector magnetogram sequence released by the HMI 
Team, starting with a description of our preliminary 
treatment of the data ($3]) followed by a step-by-step 
demonstration of the method (Q and analysis of its re- 
sults (©. We conclude with a brief summary in <J6] 

2. CALIBRATING DOPPLER SHIFTS WITH FARADAY'S 

LAW 

Consider a region of the photosphere where flux is ei- 
ther increasing or decreasing due to emergence or can- 
cellation (by any of the mechanisms discussed above), 
respectively. In an area A containing flux of one polar- 
ity, Faraday's law and Stokes' theorem relate the time 
rate of change of magnetic flux, d&/dt, through A with 
the integral of the electric field cE projected onto a closed 
curve C bounding A, 
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where n is the unit normal vector on A and xii is the 
tangent vector to C oriented to circulate in a right-hand 
sense (counterclockwise) with respect to n. This must 
apply separately to the areas A± containing the flux of 
each polarity. 

We restrict ourselves to the case where the flux through 
areas A± changes only due to emergence or cancellation, 
so flux is only added to or removed from A± along the 
PIL. Hence, the only contribution to the integral of E 
along C occurs along the PIL, 
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The time rate of change of flux is therefore equal to 
the voltage drop along the PIL. This equivalence is ex- 
act, without approximation. In real data, however, PILs 
might not be observable at a given instrument's spatial 
resolution, flux cannot be measured without uncertain- 
ties, and flows on the peripheries of A± can disperse flux 
until its density falls below a magnetograph's sensitivity. 

Whether the electric field contains a non-ideal compo- 
nent or not, the rate of ideal transport of magnetic flux 
across the PIL is given by replacing cE with — (v x B) 
in equation ([5]), 
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represents the transport of flux across the PIL by plasma 
flows. We now define a unit vector perpendicular to the 
PIL, xj_, with respect to n and xii, 

Xx = Xii X fl , (10) 

which points in the direction of the horizontal gradient 
of B n evaluated at the PIL (because the PIL is assumed 
to be a zero contour of i? n ), and cfxii = XiiefL. Then 
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where B^ the field tangent to the surface, B± is the 
component of B^ along xj_, and the final equality holds 
because B n vanishes along the PIL. 

If the magnetic evolution is ideal, then the rate of 
change of flux in each area A± will match the rate of 
ideal flux transport across the PIL, 
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where dA± denotes integration over cither polarity (not 
both). The key point here is that the PIL-integrated 
rate of transport of horizontal flux by the normal veloc- 
ity matches the rate of change of normal flux through 
A±. Unlike equation (|5|). this equivalence rests upon 
the assumption of ideality. In general, only the unsigned 
fluxes match, since both up- and downflows can lead to 
flux decrease by cancellation (Figure [TJ bottom left and 
center), and both up- and downflows can lead to flux 
appearance, by either fi-loop eme rgence or f7-l oop sub- 
mergence (Figure [TJ top row; and lAbbettl[2Q07h . Fi sure 
[3] sketches emergence of an fi-loop flux system, one phys- 
ical situation to which this formalism can be applied. 

We can apply this formalism to PILs of the LOS field 
near disk center. We require both Dopplergrams, in 
which the LOS velocity t>LOS is measured, and vector 
magnetograms, in which pixel-averaged flux densities of 
magnetic components both along the LOS, -BloSi and 
transverse to the LOS, B trs , are determined. We approx- 
imate the normal vector n used above with the unit vec- 
tor along the LOS, £, replace the exact time derivatives 
above with finite-difference approximations, and change 
integrations to summations over pixels in A± and along 
the PIL. Equation ([S]) then becomes 
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where Ax is the pixel length, and each sum — there 
are two, one over each polarity — runs over pixels in 
the neighborhood of the PIL. (We address identification 
of PILs and definition of neighborhoods near each PIL 
below.) Equation (fT"2")) then becomes 
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where the subscript vBL denotes that this expression where the sum runs over PIL pixels, and B± now refers 
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Figure 3. A depiction of ideal flux emergence across the pho- 
tosphere (gray surface). Magnetic field lines 1,2, and 3 extend 
from the solar interior (red dashed lines), across the photosphere 
(red X 's), and out into the outer solar atmosphere (red solid lines). 
Along a PIL (dashed black line) between regions of positive (white) 
and negative (black) magnetic flux, velocity (green vector) normal 
to the photosphere is driving flux emergence. At the instant de- 
picted, the field line labeled 4 is crossing the photosphere. Ad- 
ditional vertical transport of flux (containing field lines 5 and 6) 
would cause unsigned flux through the positive and negative flux 
regions to increase. An ideal electric field (blue vector) parallels 
the PIL at the emergence site. 



to the component of B trs along xj_ . Equation (|14|) then 
becomes, formally, 
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where we have multiplied equations (fTS"]) and ([TBI by At 
to deal with changes in flux instead of rates of change in 
flux, and the sum in the left half of equation (|X8|) runs 
over the near-PIL neighborhood of one LOS polarity or 
the other. 

Crucially, LOS velocities along PILs are therefore con- 
strained by evolution of nearby LOS flux: over a time 
interval At, the change in LOS flux in each polarity in 
the neighborhood of a PIL should match, within method- 
ological uncertainties, the transport of transverse mag- 
netic field by LOS velocities summed along PIL pixels. 

As mentioned in the introduction, however, a nonzero 
bias velocity vo might be present in estimated Doppler 
velocities from HMI, perhaps due to instrumental effects 
or the convective blueshift. We define the biased and 
true Doppler velocities as v^ os and vloSj respectively, 
which are related via 



J LOS 



"LOS + Vq 



(19) 



we follow the astrophysical convention that receding ve- 
locities with respect to the observer — i.e., redshifts — 
are positive. Then for v > 0, the biased Doppler shift 
is more red than the true Doppler shift. Such an offset 
in the Doppler velocity would mean that the true rate of 
flux emergence or submergence is related to the biased 
rate At>' vBL /At by 



A$, 



BL 



At 



J2 Axv 'losB±-J2 Axv o b ^ ( 2 °) 



At At ' (21) 

where the bias flux A$bias is due to the bias velocity vq, 



A$v 



v AtBL , 



(22) 



and where we define the magnetic length of the PIL as 



BL 



PIL 



(Ax)B A 



Equation (|17[) then implies 
|A$T,ns. + l = IA$' 



vBL 



(23) 



(24) 



In the absence of errors, the observed quantities A$los 
and A$' uBL enable A<&bi as to be determined, from which 
the bias velocity Vq can be found via equation (|22|) . Er- 
rors are certainly present in these quantities, however, 
and the evolution also might not be ideal. At this point, 
however, two key points should be noted: 

• First, the rest wavelength of the spectral line used 
to infer Doppler velocities is a unique, well-defined 
physical quantity. 

• Second, any bias in the rest wavelength should be 
spatially uniform over the instrument FOV in an 
individual measurement. 

These ideas imply that the set of bias fluxes {A«f>bias} and 
magnetic lengths {BL} measured over a set of PILs near 
disk center can be used to estimate vo statistically, which 
also enables quantifying uncertainties in the estimate. 

As noted in the introduction, away from disk center, 
the approximate coincidence between LOS PILs and the 
radial-field PILs where flux appears or cancels breaks 
down. While the formalism used here could potentially 
be developed further to enable analysis of PILs signifi- 
cantly away from disk center, we will restrict our analysis 
here to PILs near disk center. In Appendix [Xl we con- 
sider pathologies in applying this method to PILs away 
from disk center. 

There are several ways one might estimate the flux 
changes A$los and A$ v bl and magnetic length BL 
for each PIL, along with several ways sets of these esti- 
mates can be used to estimate vq, and thereby remove 
any bias present in measured Doppler velocities w L os- 
Accordingly, in the following sections, we demonstrate 
our approach with actual HMI observations, with sev- 
eral goals in mind: first, to make the method more clear; 
second, to show that the method can be used with real 
data; third, to show the method derives physically rea- 
sonable values; and fourth, to investigate the extent to 
which the method depends upon input parameters and 
is susceptible to errors. 

3. HMI OBSERVATIONS 

Recently, the HMI Team released a sequence of 
"cutout" vector magnetogramsQ and DopplergramtQ 
from NOAA AR 11158, from 12-16 Feb. 2011. Members 
of the HMI Team have put detailed information about 

1 ftp:/ /pail. stanford.edu/pub/HMIvector/Cutout/ 

2 ftp:/ /pail. stanford.edu/pub/xudong/stage2/IcV/ 



PIL 



PIL 



7 



this dataset online^, and are preparing a paper describ- 
ing production of this dataset. This active region was 
the source of an X2.2 flare on on 2011/02/15, starting in 
GOES at 01:44, ending at 02:06, and peaking at 01:56. 
Figure 2] shows the LOS magnetic field in a subregion of 
the data array near 10:00UT on 2011/02/14. 

In this section, we briefly describe procedures that we 
undertook to prepare this data prior to applying our cal- 
ibration procedure^ 

After downloading the FITS data files, several pro- 
cessing steps were required, which we describe here. 
Because the reacLsdo .pro IDL procedure did not (at 
the time of this writing) properly handle pixels set to 
the BLANK value in the FITS headers, we used the 
f itsio_read_image .pro procedural to read in the data. 

3.1. Cropping 

To reduce the full dataset to a more manageable size, 
we focus on a subset of the full five-day sequence. To 
study photosphcric magnetic evolution prior to the X 
flare, and to baseline this evolution against post-flare 
evolution, we retain about 72 hours of data, from mid- 
night at the start of the 13th until midnight at the 
end of the 15th, inclusive of endpoints. Given the 12- 
minute cadence, the resulting time series consists of 
361 time steps. The active region was at S19E11 at 
00:30UT on 2011/02/13 and S21W27 at 00:30UT on 
2011/02/16. Pixels lacking data near the edges of the 
cutout FOV, from artifacts of the cutout process, were 
cropped: columns [0 - 24] and rows [0 - 5] were removed 
for all steps. 

3.2. Removing SDO Motion and Solar Rotation 

Next, Doppler velocities were corrected for spacecraft 
motion. Due to SDO's geosynchronous orbit, its velocity 
along the radial, Sun-observer line can be large, of order 
± 5 km s^ 1 . Also, because the Earth is orbiting westward 
about the Sun, and SDO orbits Earth, there is a signif- 
icant projection of SDO's westward motion onto lines of 
sight to many pixels. In addition, there is a nonzero 
component of SDO's northward velocity onto the LOS. 
While the radial component is larger than the W and N 
components, the latter can be significant — a few tens 
of m s _1 or more. Accordingly, the spacecraft velocity 
was projected onto the LOS to each pixel — in the small- 
angle approximation — and subtracted off the measured 
Doppler velocities. Then Stonyhurst latitudes and lon- 
gitudes of each pixel (accounting for the solar B and P 
angles) were computed to correct Doppler velocities for 
solar rotation , using the "magn etic / 2-day lag" rotation 
rate found by iSnodgrasi (|1983fl , for which 



Wrot = A + B sin(0) 2 + C sin(6>) 4 
A = 2.902 
B = -0.464 
C = -0.328 , 



(25) 
(26) 
(27) 
(28) 



3 http:/ /jsoc. stanford.edu/jsocwiki/VectorPaper 

4 More detailed notes de scribing our procedures a re online at 
http : / /solarmuri . ssl. berkeley.edu /~ welsch / public /manuscripts / 
Doppler _cali b / hmi_data_notes_current . txt 



5 Available online at http://www.mps.mpg.de/projects/seismo/ 
GDC_USE/using_drms.html, along witn a required, compiled 
shared-object file, fitsio.so. 



where w rot is the photospheric rotation rate, 9 is lati- 
tude, and A,B, and C are in microrad s _1 . Accurate 
correction for rotation prior to application of our cali- 
bration procedure is not necessary; but the bias velocity 
estimates will include (and could be dominated by) the 
Doppler shift from rotation. 

3.3. Removing Azimuth Reversals 

Accurate estimation of electric fields by PTD (or, re- 
ally, any method) requires that changes B between mag- 
netograms arise from physical processes on the Sun, as 
opposed to artifacts of the measurement process. Hence, 
any spurious changes in the measured field, where they 
can be identified, should be removed. 

We noticed that in successive images of the transverse- 
field azimuths, the inferred azimuths in some regions 
flip by nearly 180 degrees from one frame to the next, 
as would be expected from errors in resolving the 180- 
ambiguity. The seven successive frames (12 min. apart) 
in the left column of Figure [5] show an example of this 
effect, with time increasing downward. The area shown 
corresponds to that in the solid white box in Figure 21 
and the initial frame corresponds to that in Figure 2J 
Azimuths range over [0,360] degrees. The next column 
shows maps of the angular difference (the interior angle) 
between the current and previous azimuths. Differences 
in azimuths range from [-180,180] degrees (note the ro- 
tated color wheel in the top panel) . The black patches in 
rows 4 - 6 of this column correspond to a region in which 
azimuths flip by nearly 180° in a region from one frame 
to the next. The checkerboard pattern in the same region 
in rows 2 - 3 of this column implies a spatially alternat- 
ing pattern of 180° azimuth flips, and probably arises 
from lack of convergence in the simul ated annea l ing a l- 
gorithm used to infer the azimuths ([Leka et al J 120091 1 . 
Also, the red-and-blue, speckled regions toward the left 
side of each frame in this column exhibit rapid spatial 
variations in azimuth changes. These speckled regions 
correspond to areas with relatively weak field strength 
(B < 200 G), where azimuth determinations could be 
more problematic than in strong-field regions. 

Interpreting frame-to- frame azimuth flips over finite re- 
gions as spurious, we seek to identify and remove these 
artifacts. Our goal is to automate detection of such flips, 
and our basic approach is to identify suspicious changes 
in azimuth, which we envision as "top hats" (or inverted 
top hats) in the running differences of azimuths in in- 
dividual pixels: large, positive (or negative) jumps in 
azimuth for one time step, followed by reversals to the 
pre-jump level at the next step. The actual patching is 
simpler: we add 180 degrees to the azimuth, and output 
the flipped data modulo 360 degrees. 

We tried more than one approach to detecting and cor- 
recting spurious azimuth changes before finding one that 
we think works well. To save other researchers from re- 
peating our efforts, we now describe, in detail, both our 
failed approach and the approach we found worked best. 

In our initial, unnsuccessful approach, we simply iden- 
tified all pixels with unsigned changes in azimuth from 
frame (i— 1) to frame i above a given threshold (we tried, 
e.g., 120, 135, 150, and 170 degrees), and flipping these 
in frame i. The updated frame i was then used as a refer- 
ence for finding large changes in frame (i + 1). Using this 
approach, once a pixel is flipped, its new state tends to 
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Figure 4. A grayscale image of the LOS magnetic field in AR 11158 near 10.00UT on 2011/02/14 in subregion of the data released by 
the HMI Team. White/black represent positive/negative LOS flux, respectively. Image saturation is set at ±750 G. The solid white box 
outlines a subregion of the FOV shown in Figure\5\ and the white/black dashed box outlines the subregion shown in Figure^ 



be propagated forward in time. But legitimate magnetic 
evolution can cause this state to be incorrect, leading to 
spatially isolated pixels that differ from their neighbors. 
This contradicts a central tenet of the original ambigu- 
ity resolution procedure: currents should be minimized. 
Further, the number and spatial arrangement of prob- 
lematic pixels appears to depend strongly on threshold. 

The best approach we found is to flip azimuths only in 
pixels where both: (i) the jump in azimuth was larger 
than 120° between frames; and (ii) the jump in az- 
imuth increased the time-averaged, unsigned changes in 
azimuth over a the set of frames within ±N steps of a 
given frame. Use of a range of images on either side 
of the i-th image is necessary to capture the essential 
character of a top hat: a large change, followed by a re- 
versal. (A large azimuth change from one frame to the 
next, by itself, might be legitimate — e.g., by horizon- 
tal advection of a region with a spatial transverse field 
reversal.) We determined suspicious flips by comparing 
means of two arrays of 2N azimuth differences: (1) un- 
signed differences between the azimuths at frame i and 
azimuths in ±N neighboring frames in time, [i— N : i — 1] 
and [i + 1 : i + N); and (2) unsigned differences be- 
tween the flipped azimuths (hypothetical data, with all 
pixels' azimuths flipped) at frame i with actual azimuths 
in ±N neighboring frames in time, [i — N : i — 1] and 
[i + 1 : i + N]. Pixels for which flipping would decrease 
the mean, unsigned, frame-to-frame angular differences 
are then flipped. 

We call frames in which azimuths have been flipped 
filtered. Filtered azimuths are used for past frames, 
[i — N : 2 — 1], while unfiltcrcd azimuths are used for 
future frames [i + 1 : i + N] , so this approach is not time- 
symmetric. This approach should be able to deal with 



top hats that are N steps wide, but not wider. We tried 
both N = 4 and N = 2. Using N = 4 results in 2-5% 
fewer flips in a given frame, so is slightly more restrictive. 
We did not, however, see much difference between results 
from N = 2 and N — 4 in strong-field regions, so most 
differences are probably in weak-field regions where in- 
ference of the direction of transverse fields is less reliable 
anyway. 

The third column from left in Figure [5] shows filtered 
data from the same seven successive frames as in the left- 
most column, with suspicious changes in azimuths iden- 
tified by our approach flipped by 180°. The right-most 
column shows maps of the angular difference between the 
current and previous azimuths (note the rotated color 
wheel in the top panel) from the third column. The 
checkerboard and black patches visible in the second col- 
umn are not present in this column. Also, changes in 
azimuths in the red-and-blue, speckled regions on the 
left sides of each frame are not as large in this column 
as in the second column. This suggests our approach 
decreases fluctuations of azimuths in weak-field regions, 
too. 

4. DEMONSTRATION OF ELECTROMAGNETIC 
CALIBRATION 

We now seek to quantify any bias velocity in the 
Doppler velocity measurements. Since this will require 
associating pixel values in successive frames (e.g., in 
equation [T5|) . we first co-aligned the plane-of-sky (POS) 
Dopplergram, field strength, and inclination arrays. We 
used the £?los arrays as the reference obser vations, since 
struc tures in -Blos are long-lived (e.g., iWelsch et al.l 
l2012f h Whole- frame shifts between each -Blos array were 
determined to sub-pixel accuracy using Fourier cross- 
correlation, and the data arrays were shifted accordingly 
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Pixel 

Figure 5. Wc applied a filtering procedure to remove back-and- 
forth dips in azimuth, which are likely due to frame-to-frame incon- 
sistencies in ambiguity resolution. From left to right, 1st column: 
100 X 100 pixel maps of transvcrsc-Rcld azimuths (color wheel in 
top panel ) from the solid white box in Figure [4] from seven suc- 
cessive vector magnetograms (separated by 12 minutes), starting 
at 10:00UT on 2011/02/14. Time increases down. 2nd column: 
Difference in azimuth between current and previous maps (note 
different color wheel in top panel). Initial difference map is arbi- 
trarily set to zero. 3rd column: Filtered versions of the azimuth 
maps in the 1st column, in which continuity in azimuths has been 
imposed (color wheel in top panel). 4th column: Difference in az- 
imuth between current and previous filtered maps (color wheel in 
top panel). Initial difference map is arbitrarily set to zero. Back- 
and-forth flips in many pixels have been removed. 

via Fourier interpolation^ Instead of wrapping data 
across image edges, due to the assumed periodicity in 
the Fourier method, data shifted out of the image field 
of view were zeroed out. Interpolated field inclinations 
outside [0,180] were capped at these values. We then used 
the co-aligned inclinations and field strengths to derive 
co-aligned LOS and transverse fields. 
Our calibration method requires several additional 

6 See http://solarmuri.ssl.berkeley.edu/~welsch/public/ soft- 
ware/shift_data. pro- 



tasks: identifying PILs of the LOS field; quantifying 
changes in LOS magnetic flux near those PILs; summing 
transverse fields B trs arid Doppler velocities ^los a l° n g 
identified PILs; and statistically estimating any offset vq 
in these Doppler velocities. 

4.1. Identification of PILs 

The first step in our approach is to use an automated 
algorithm to identify PILs in LOS magnetograms. Our 
procedure creates masks of each polarity for pixels with 
|-Blos| above a field strength threshold of -Bthr, dilates 
each mask by one pixel, and finds all areas where the 
dilate d masks overlap (see lSchriiverll2007l : lWelsch and Lil 
2008). The threshold field is the only input parameter in 
PIL identification procedure; here, we use Bthr = 50 Mx 
cm~ 2 , ensuring that we identify relatively strong fields 
in close proximity. This approach identifies structures 
1-3 pixels wide, which we erode into single-pixel-width 
lines with IDL's morph_thin.pro procedure. We define 
the resulting pixels to be PIL pixels. In Figure H3 wc 
show identified PIL pixels, color-coded by (possibly bi- 
ased) Doppler velocity i^Qg, m a c l° se ~ u P view of the 
LOS magnetic field at the time step when AR 11158 was 
closest to disk center. 

PILs & Uncorrected LOS Velocities, 20 1 1 .02. 1 4_02:59:54 
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Figure 6. Grayscale image shows a zoomed-in view of £?LOS> 
with PIL pixels superimposed and color-coded by (possibly bi- 
ased) Doppler velocity v' l q S . Grayscale saturation is set to 100 
Mx cm~ 2 , and the velocity saturation is set to ± 250 m sec -1 . 
JVote the predominance of redhsifted PILs. 

In this sub-region of the full active region, iVpiL = 
98 PILs were identified. Most of the Doppler signals 
along PILs appear redshifted. These are almost cer- 
tainly pseudo-redshifts, arising from the Doppler zero- 
point being defined to be blue- ward of the true rest wave- 
length for photospheric plasma with no LOS velocity. 
The velocity zero point for HMI is derived from the me- 
dian of Doppler velocities within 90% of the solar radius 
(Y. Liu & S. Couvidat, private communication, 2011). 
This implies contributions from quiet-Sun (and therefore 
strongly convecting) plasma dominate the median, and 
probably bias the zero-point. 

4.2. Flux Changes Near PILs 

Having identified PILs at a given time t , our sec- 
ond step is to estimate the set of changes of LOS flux 
{A«f>LOs} near all PILs over the At = 24 minutes be- 
tween the two LOS magnetograms observed ±12 min- 
utes from the magnetogram at to. To do so, for each 
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PIL we create a binary map of that PIL's pixels, then 
apply IDL's dilate .pro procedure with a (d x d) matrix 
of l's to the binary map to define that PIL's "neigh- 
borhood mask." The dilation parameter d is the only 
other free parameter in our method. While granular 
flows have speeds on the order of ~ 1 km sec -1 , such 
flows are short-lived, with lifetimes of ~ 5 min. Over 
At = 24 min, therefore, displacements of fluxes tend to 
be consistent with lower averaged speeds, o n the order of 
~ 0.5 km sec -1 or less (see, e.g., Fig. 13 of iWelsch et al.l 
l20ll . Given the HMI pixel size Ax ~ 0.5" ~ 360 km 
near disk center, magnetic flux should not move by more 
than about two pixels over At = 24 min. Here, we use 
d = 7 — that is, dilation by ±3 pixels in all directions 
— although we have explored other choices (see below). 
Next, we multiply the two co-aligned LOS magnetograms 
from ±12 minutes by this neighborhood mask and sum 
the unsigned flux in each product, then difference these 
summed fluxes. Because this counts LOS flux changes in 
both polarities, we divide the result by two to estimate 
A$los over At — 24 minutes. 

It should be noted that converging or diverging hor- 
izontal flows unrelated to emergence/ submergence at 
each PIL can transport "background" flux into the neigh- 
borhood of a PIL, and near-PIL flux out of this neigh- 
borhood, in some cases dominating our estimates of the 
change in LOS flux due to emergence or submergence 
along the PIL. Hence, the computed rate of change in 
LOS flux along an individual PIL is subject to large er- 
rors. Identifying collections of like- polarity pixels as "fea- 
tures" and tracking their evolution (|DeForest et al.ll2"o"07t 
IWelsch et al.ll20lili might be one way to better account 
for fluctuati ons in LOS flux due to flux concentration 
/ dispersal (|Lamb et al.l 120081 ). Assuming the contribu- 
tions from converging and diverging flows cancel in the 
aggregate over many PILs, the statistical properties of 
changes in LOS fluxes along a set of PILs are more ro- 
bust to such errors, and therefore more useful for our 
purposes. 

Our third step is to compute, for each PIL, the corre- 
sponding expected flux change from the (assumed biased) 
Doppler signal, A<fr' vBL , and magnetic length, BL. Since 
the expressions derived for these quantities (equations 
IT8l and |23|) in S}2] refer to the component of the trans- 
verse magnetic field perpendicular to the PIL, we inves- 
tigated the angle of the transverse magnetic field along 
PILs in our dataset. Figure [7] illustrates the geometry of 
the transverse field along several PILs. The upper-left 
panel shows unit vectors perpendicular to the PIL, 
determined from the gradient of -Blos evaluated at PIL 
pixels. The upper-right panel shows unit vectors b tIS 
along B trs , evaluated at PIL pixels. Note predominance 
of unit vectors b trs pointing toward the bottom of the 
FOV: this coherence might be spurious, and could arise if 
the minimization of currents by the ambiguity resolution 
algorithm introduces artificial, large-scale correlations in 
azimuths. Color coding in the bottom-left panel shows 
the cosine of the angle x between these unit vectors — 
their dot product, xj_ • b trs — with red positive, and 
the saturation level set to ±0.866. (The sign of the an- 
gle x is defined such that the angle between xj_ and X|| 
is +90°; see equation ITUl') Fields along PILs evidently 
exhibit components in both the "normal" (positive-to- 



negative; b lue) and "inverse" (negatiy e-to-positive; red) 
directions (|Martens and Zwaanl l200ll ) perpendicular to 
the PIL. For each PIL, the mean (PIL-averaged) signed 
and unsigned cos(x) and unsigned sin(x) can be com- 
puted; the bottom-right panel shows the distributions of 
these averages for all PILs in our dataset. The peaks at - 
1 and 1 for the signed cos(x) and around 0.5 for | sin(x)| 
imply that the average field of many PILs: (i) is spa- 
tially coherent; and (ii) tends to point primarily across 
(not along) the PIL. 

As with changes in LOS flux near PILs — which can 
arise from converging or diverging flows unrelated to 
emergence / submergence — some component of Doppler 
velocities along identified PILs might be unrelated to 
emergence / submergence. Slight inaccuracies in PIL 
identification, for instance, could lead us to improperly 
include flows along the magnetic field (siphon flows) in 
our flux transport rates. Also, filling factors — due to 
scattered light within the HMI instrument, or imaged 
emission from unresolved, unmagnetized plasma, or both 
— could bias the estimated flux transport rate in a given 
pixel (and contributions from either would plausibly tend 
to be blue-shifted). Hence, as with changes in LOS flux 
near PILs, the inferred rate of transport of transverse 
flux along an individual PIL is subject to large errors. 
Again, the statistical properties of flux transport rates 
along the aggregated collection PILs should be more ro- 
bust to error, and therefore more useful for quantifying 
the pseudo-redshift. 

We have also investigated making the simplifying re- 
placement B± —> B trs when computing AQ v bl and BL. 
Physically, one can justify this approximation by noting 
that the footpoints of fields with components tangent to 
the PIL must still thread the photosphere somewhere, if 
not directly across the PIL, as would be the case for no 
magnetic component tangent to the PIL. We will show 
results computed both ways, using either B± or B tTS at 
each PIL pixel, but unless otherwise stated, results shown 
below are derived using i?trs- As will be seen, this did 
not drastically change the estimates of bias velocities vq, 
although quantities derived using B± were substantially 
noisier. We also note that using J5 trs makes the bias es- 
timation technique independent of ambiguity resolution. 

4.3. Estimation of Bias Velocities 

Finally, we can estimate any bias velocity present. We 
first difference A<I>los — AQ' vBL on each PIL to com- 
pute the set of bias fluxes {A$bias} from all PILs. We 
must then esti mate the coefficient vq from the ratio of 
{A$ bias } to {AtBL}. 

In the data, errors are present in both measured quan- 
tities. Before estimates of typical uncertainties were pub- 
lished by the HMI team, we adopted uniform uncertain- 
ties of 20 m/s for the Doppler velocities, 25 Mx cm~ 2 
for Blos j and 90 Mx cm~ 2 for B tls . This value for 
the Doppler uncertainty is consistent with expected near - 
disk-center noise levels reported by ISchou et al. (2012). 
The value for noise in -Blos is significantly larger than 
the HMI Team's since-published estimate the uncertainty 
of ~ 6 Mx cm~ 2 in 720-second data (e.g., Liu, Scherrer, 
et al. 2012), while the value for B trs is probably closer to 
the HMI Team's estimate of noise in t he transverse com - 
ponent on the order of 10 2 Mx cm~ 2 (|Sun et al.ll2012D . 
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Figure 7. Top left: Unit vectors normal x± to PILs, along the gradient in -Blos along each PIL in a small subregion of the full FOV. 
Grayscale background shows Blos> with saturation at ±100 Mx cm -2 . Top right: Unit vectors btrs parallel to Btrs along PILs in the same 
FOV. Bottom left: Color-coded dot product of x± ■ btrs, with red positive ("inverse" polarity) and saturation set to 0.866, corresponding 
to 30° . The component of Btrs perpendicular to each PIL is homogeneous along many PILs. Bottom right: Distributions of PIL-averaged 
cosx (black), |sinx| (red), and |cosx| (blue), where \ is the angle between and btrs- Fields tend to point across (not along) most 
PILs. 



We estimate vq and its uncertainty by analyzing the 
the set of ratios {A$bias/(Ai-BL)} over a subset of iden- 
tified PILs. Based upon both our prior knowledge that 
the convective blueshift biases the median-derived esti- 
mate of the rest wavelength in HMI Doppler data blue- 
ward, and the predominance of redshifted PILs in Figure 
[51 we expect that the biased PIL velocities are more red 
than the true velocities. Recalling that we use the astro- 
physical convention that redshifts correspond to positive 
velocities with respect to the observer, the bias velocity 
vo should be positive. 

Not all bias fluxes, however, arc consistent with a posi- 
tive bias velocity. This is not obvious from equation (|24[) . 
since it deals with absolute values. Consequently, we now 
consider the different possibilities for the relative sizes of 
A<I>los, A&' vBL and A$bias- A positive vo is consistent 
with PILs that obey either 



or 



A<P' vBL > |A$ LOS | 



|A$ LOS | < A$' uBi <0 



(29) 



(30) 



since the correction to A$' 
For PILs that obey either 



BL 



is -A$L M = -v AtBL. 



|A$los| 



> a& vBL > 



or 



a$: 



vBL 



< 



|A$ 



LOS 



< 



(31) 



(32) 



however, vq would have to be negative to improve agree- 
ment between the flux changes. 

What fraction of PILs are consistent with a positive 
bias velocity vq? Of the ~50,000 PILs identified in the 
full dataset (with B thr = 50 Mx cm~ 2 ), 73.8% and 5.7% 
were consistent with equation (|2"9"j) and ([50)1 . respectively, 
and 13.9% and 6.6% obeyed equations (|3"T|) and 
respectively. Consequently, ^80% of PILs are consistent 
with a positive bias velocity. 

In Figure [51 we show the distribution of bias fluxes for 
all PILs at all time steps. (In any single time step, there 
are ~ 140 PILs in our FOV, too few to form a continuous 
distribution.) The distribution is strongly skewed toward 
positive bias fluxes. We consider PILs matching either 
(|3Tj) or (|32|) to be pathological, which we attribute to 
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"noise," primarily systematic errors in our estimates of 
the fluxes. (Errors from our assumed uncertainties in the 
underlying magnetic and Doppler data would be 3 x 10 16 
Mx per pixel or less.) 

This distribution of bias fluxes can be used to quan- 
tify uncertainties the in bias flux. Mirroring the negative 
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Figure 8. Dashed line shows the total distribution of bias fluxes 
for all PILs at all time steps, which is clearly skewed toward the 
positive. Mirroring the negative bias flux values across zero and 
fitting the result with a Gaussian (thin solid) enables an empirical 
estimate of noise in the bias flux from the fitted width, 3.75 X 10 17 
Mx. Subtracting the fit from the total distribution of bias fluxes 
leaves the residual distribution (thick solid) of bias fluxes above 
the noise level. 



bias flux values across zero and fitting the result with a 
Gaussian enables an empirical estimate of noise in the 
bias flux from the fitted width, 3.75 x 10 17 Mx. Fol- 
lowing the same procedure but using B±_ instead of B trs 
to compute A$' vBL results in a fitted width of 5 x 10 17 
Mx, reflecting greater uncertainty in estimates of this 
flux. These define thresholds in the bias flux: estimates 
of bias fluxes smaller than these levels fall within our 
(systematic) uncertainties. 

Subtracting the fit from the total distribution of bias 
fluxes leaves the residual distribution of bias fluxes above 
our uncertainty level. At each step, we computed sepa- 
rate estimates of the bias velocity from both the full set 
of PILs and just the subset of PILs with bias fluxes more 
positive than these estimates of the threshold bias fluxes. 
Not surprisingly, estimates of vq using only PILs with 
bias fluxes above these thresholds — which we refer to 
as the "high-bias" subset of PILs — were substantially 
higher (as shown below) than estimates derived using all 
PILs. We view use of all PILs as more conservative, so 
results shown below are derived from all PILs, unless 
explicitly stated otherwise. 

5. RESULTS & ANALYSIS 



In Figure UH we show sorted values of the set {vq} = 
{A$bias/(Ai-B£)} from all PILs, scaled to units of m/s, 
derived from three successive magnetogram triplets when 
AR 11158 was close to the central meridian. It is clear 
that most values fall well above zero (the lower, dotted 
line in each plot; recall that we define rcdshifts as posi- 
tive). This is entirely consistent with the predominance 
of redshifts in Figure [6] The means and standard er- 
rors in the estimates are 101 ± 17 m/s, 126 ± 18 m/s, 
and 138 ± 16 m/s, respectively; the means are plotted as 
horizontal solid lines. Medians are the dashed horizontal 
lines in each plot (indistinguishable from the mean in the 
bottom plot). We use the standard error in the estimate 
(standard deviation divided by %/^pil) as a measure of 
uncertainty since each estimate is a measurement of the 
same number, vq, although we do not know if the errors 
are Gaussian. Only ~ 45% of error bars are consistent 
with the average vq in each plot. This suggests either: 
(1) that our error estimates are too low, perhaps due to 
our neglect of systematic errors (e.g., in definitions of 
PIL and their neighborhoods, and in quantifying mag- 
netic fields along PILs); or (2) that electric fields along 
PILs are in some cases inconsistent with the ideal electric 
field assumed in equation 
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Figure 9. Sorted values of the set of bias velocities, {vq} = 
{ A$ bias /( AtBL)} from all PtLs, (in m/s) and estimated errors 
for each PtL (with the astrophysical convention of positive red- 
shift), from three successive (and partially overlapping) magneto- 
gram triplets. The lower thin dotted line shows the (uncalibrated) 
zero Doppler velocity. Solid and dashed horizontal lines show av- 
erage and median bias velocities. The magnitudes of the estimated 
bias velocities are 101 ± 17 m/s, 126 ± 18 m/s, and 138 ± 16 m/s 
average velocity, suggesting our error estimates are too low, or the 
ideal assumption is invalid. Nonetheless, our results are consis- 
tent with the predominance of redshifts in Figure [6) and indicate 
pseudo-rcdshifts are present. 
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In the top panel of Figure [TU] we plot estimates of vq 
versus time, for different approaches — using either B trs 
or B±, and either all PILs or only high-bias PILs - 
along with different choices for the two parameters in 
our method, the threshold LOS field strength B t \ lr (25 
or 50 Mx cm -2 ) that defines significance in identifying 
PILs, and the dilation parameter d (5 or 7) that defines 
the neighborhood near each PIL over which changes in 
LOS flux with time are computed. Estimates made from 
the median vo from all PILs, with B t h r = 50 Mx cm~ 2 
and d = 5 using either B trs or B± (red, solid and dashed, 
respectively) are very similar. (We use medians instead 
of averages here because some estimates in the full set 
of PILs made using B± are pathologically noisy.) The 
blue lines are from averaging only high-bias PILs, using 
either B trs or B± (solid and dashed, respectively), and 
also agree closely. Using B trs , increasing d to 7, with 
-Bthr = 50 Mx cm" 2 , gives the thick, black dashed line 
(which closely follows the red lines); decreasing B t h r to 
25 Mx cm -2 , with d = 5, gives the blue dash marks 
(which are also close to the red lines). Together, these 
suggest that results are not strongly affected by varying 
d slightly, or using B trs versus B±, although estimates 
of vq from high-bias-flux PILs are systematically higher. 
We note that the less restrictive threshold flux density, 
Bthr = 25 Mx cm~ 2 , resulted in somewhat lower esti- 
mated bias velocities, especially toward the end of the 
run. This might be because PILs in weaker-field regions, 
which should be more strongly convecting, were included 
when the lower threshold was used. 

The HMI Team's estimate of the 1 — a noise in Blos 
is about 6 Mx cm -2 , so the 50 Mx cm -2 threshold corre- 
sponds to about 8 standard deviations. With Bthr = 50 
Mx cm -2 and d = 5, the standard deviation between the 
i-th bias estimate made using B trs , and a 5-step (1-hr.) 
running average was 18 m/s, near the noise level of the 
Dopplcrgrams. Standard errors in the mean (SEM) were 
computed using uncertainties assumed above; the aver- 
age SEM for this series was 22 m s _1 . The consistency 
of the vq estimates on hour-long time scales is evidence 
of robustness in our estimates. Variation of vo on time 
scales longer than few hours is evidence that calibration 
of the Doppler shifts in time is necessary, as the estimated 
bias velocity is not constant. For comparison with the 
phase of SDO's orbit, we also plot the radial component 
of spacecraft velocity, rescaled and shifted in the vertical 
direction (but not in time) , in the sinusoidal, black, solid 
line. 

The solid, black, jagged line shows an empirical esti- 
mate of vq: the median of Doppler velocities along all 
PIL pixels at each step, identified with the 50 Mx cm~ 2 
threshold. This estimate is relatively simple to deter- 
mine: one only needs to identify PILs of the LOS field, 
then take the median Doppler velocity on all PIL pixels. 
This approach, however, ignores the flux-matching con- 
straint in equation (fT8")) . and therefore could be biased 
by a strong episode of flux emergence (or submergence) . 

A trend in each curve with longitude can be noted, of 
about 5 m s _1 deg -1 for the curves shown. In terms of 
an error in our assumed rotation rate, this corresponds to 
~ 0.45/xrad s _1 , a significant fraction of the A coefficient 
in equation (|25p . We note that accurate compensation 
for the rotation rate prior to applying our method is not 



required: it can also be applied even if the rotational 
Doppler shift has not been removed, although the de- 
rived bias velocity will then include both rotational and 
convective shifts. 

All of our estimates of the bias velocity vq are positive, 
although there is a significant spread in the estimates, 
which range from 50 - 500 m s _1 , with variations in 
both time (and therefore longitude) and from method 
and parameter selection. We will compare the physical 
implications of the differing estimates below. 

Could an episode of strong flux emergence somehow 
influence our estimate of the bias velocity? In princi- 
ple, the flux-matching constraint in equation (|18[) should 
make our approach insensitive to the rate of flux increase 
or decrease: changes in LOS flux and the transport of flux 
along the LOS should be consistent, regardless of emer- 
gence / submergence rates. Nonetheless, could system- 
atic errors in our approach make our calibration method 
susceptible to bias during episodes of strong emergence? 
We investigate this possibility in Figure [TUJs bottom 
panel, which shows the total unsigned radial flux in pix- 
els with absolute radial flux density above 50 Mx cm~ 2 , 
and the raw and smoothed (using a five-step boxcar av- 
erage) finite-difference time rate of change in that flux, 
which is positive when flux is emerging. The smoothed, 
de-meaned, and rescaled bias velocity derived from all 
PILs with B t hr = 50 Mx cm -2 and d = 5 is also plotted 
(red line). Although variations in the rate of change of 
flux do not appear to strongly influence the bias velocity 
estimates, we find a weak correlation (~ 0.2) between 
these time series. A scatter plot and linear fit, in Fig- 
ure [TTJ confirms that the dependence is weak. The fitted 
slope is -9 m s" 1 (10 17 Mx s -1 ) -1 , implying that typi- 
cal fluctuations of AQ/At, on the order of 2 x 10 17 Mx, 
should influence the estimated bias velocity by about 20 
m s , the noise level of the Dopplcrgrams themselves. 
As noted above, this correlation might arise indirectly 
from some systematic aspect of our method. An alter- 
native explanation for this correlation is the known in- 
strumental correlation between the spacecraft's periodic 
Doppler motion and the observed periodicities in HMI's 
LOS field strengths (Liu, Hoeksema, et al. 2012b). . 

In Figure [T^l we show Doppler shifts along identified 
PILs in two magnetograms, recorded about eight hours 
apart (left and right columns), both uncorrected (top 
row) and corrected by subtraction of different estimates 
of vq made in different ways: the average vq from all PILs 
(second row); the median Doppler velocity of all PIL 
pixels (third row); and the high-bias estimate of vq, using 
only PILs with bias fluxes above 3.75 x 10 17 Mx (fourth 
row). While redshifts predominate along PILs in the top 
row, Doppler shifts along PILs in the bottom rows are 
more evenly balanced between red- and blueshifts. 

Our estimates of the bias velocity fall in t he ra nge of 
a few hundred m s _1 . lAsplund and Collet) (|2003l ) used 
3D, radiative MHD simulations of magnetoconvection to 
study line profiles for several spectral lines of iron, includ- 
ing some similar to HMI's Fe I 6173 A line in wavelength. 
They find convective blue shifts similar in magnitude to 
our estimates of the pseudo-redshift (around 300 - 500 
m/s), with uncertainties (50 - 100 m/s), similar to ours 
for vq. Only vp estimates from high - bias P ILs fall in the 
range found bv lAsplund and Colletl (|2003l ). but it should 
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Figure 10. Top: Values for the bias velocity vq for different approaches and parameters used to determine changes in flux near PILs. 
The median of all PILs, using Btrs, dilation parameter d = 5 pixels, and LOS field strength significance threshold -Bthr = 50 Mx cm -2 
gives the solid red line. Error bars, omitted for clarity, are on the order of the fluctuations between successive estimates. With the same 
parameters, but using Bj_, the dashed red line results. The solid and dashed blue lines correspond to estimates using the same approaches 
and parameters as each red line, respectively, but from averaging over only high-bias PILs (see text). The dashed, thick black line (which 
coincides closely with the red lines) shows the bias velocity using Bt TB , but d = 7, keeping B thv = 50 Mx cm -2 . The blue dash marks show 
the bias velocity using d = 5, with B t hr = 25 Mx cm -2 . The solid black jagged line shows an empirical estimate of vo: the median of 
Dop pler velocity along all PIL pixels identified with the 50 Mx cm -2 threshold, which ignores the flux-emergence constraint in equation 
H8\) . The smooth, sinusoidal solid line shows l/10th of the radial component of the spacecraft Doppler velocity, offset by +250 m s~ 1 
to fit on the plot. Bottom: Total unsigned radial flux (blue) and the raw and smoothed finite-difference time rate of change in that flux 
(thin gray and thick black, respectively), which is positive when new flux is emerging. Variations in the rate of change of flux are weakly 
correlated with the smoothed, de-meaned, and scaled bias velocity estimate (red) derived using all PILs, with d = 5 and -Bthr = 50 Mx. 



be noted that their estimate of the convective blueshift 
is with respect to the true rest frame, but our estimate 
of the bias velocity is with respect to the HMI Team's 
estimate of the rest wavelength, derived from the median 
of HMFs measured Doppler velocities within 90% of the 
solar radius. 

As an aside, we remark that we found Doppler struc- 
tures along PILs to persist from one HMI vector magne- 
togram to the next, i.e., with a lifetime of at least 720 s. 



The lifetimes of patterns of upflows and downflows along 
PILs have not yet been studied yet, but bear investiga- 
tion. Autocorrelation of Dopplcrgrams mi ght be useful 
for thi s purpose, similar to methods used bv lWelsch et al.l 
(|2012l) . 

5.1. Physical Processes Along PILs 

To investigate which of our estimates for vq is most 
reasonable, we now consider the physical implications of 
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Figure 11. A scatter plot of the bias velocity as a function of 
the rate of change of magnetic flux, which are weakly correlated. 
A linear fit (dashed line) with slope 9 m s~ 1 (10 17 Mx s _1 J _1 is 
overplotted. 

our results in terms of the upward and downward trans- 
port of flux across the photosphere. Recall that processes 
that increase photospheric flux are f2-loop emergence and 
[/-loop submergence, while processes that remove flux 
are fi-loop submergence and [/-loop emergence. With- 
out correcting Doppler velocities, essentially all increases 
and decreases in flux are attributed to the two submer- 
gence processes (from [/-loops and ft- loops, respectively), 
since nearly all PILs show only redshifts. Adjusting the 
observed Doppler shifts to compensate for the convec- 
tive blueshift, however, leads to a different apportion- 
ment between the four possible process, depending upon 
the applied bias velocity: PILs that, on average, show 
upward transport of transverse flux and increases (de- 
creases) in LOS flux are ascribed to f2-loop ([/-loop) 
emergence, while PILs that, on average, show downward 
transport of transverse flux and decreases (increases) in 
LOS flux are ascribed to fi-loop ([/-loop) submergence. 

In Figure I13[ we show the frequencies of these pro- 
cesses, as functions of the number of pixels along all 
identified PILs at all time steps, for three possible cor- 
rections: the average estimates of vq, derived using B tIS 
with £?tiir = 50 Mx cm -2 and d — 5 (top); the median of 
Doppler velocities on pixels in identified PILs (middle); 
and the average high-bias estimates of vq, again using 
-Btrs with £?thr = 50 Mx cm -2 and d = 5 (bottom). In 
all panels, the black curve shows the distribution of sizes 
of all PILs. Processes involving blueshifts (f2-loop and 
[/-loop emergence) are plotted in blue, while processes 
involving redshifts (f2-loop and U -loop submergence) are 
plotted in red. Processes that remove flux via cancel- 
lation ([/-loop emergence and f2-loop submergence) are 
plotted with dashed lines. Even with the average-^o cor- 
rection (top panel), submergence processes still domi- 
nate, contradicting the expectation that the observed in- 
crease in unsigned radial flux of ~ 2 x 10 22 Mx over the 
data sequence occurred via emergence from the interior. 



Removing the median Doppler velocity over all pixels in 
identified PILs at each time step gives a more reason- 
able distribution of processes (middle panel): emergence 
and submergence processes are more closely balanced. 
But the requirement of a net flux increase implies that 
emerging processes should predominate, and this expec- 
tation is only consistent with the high-bias-i>o correction 
(bottom panel). 

In principle, observations of changes in total unsigned 
radial flux in the FOV should be consistent, separately, 
with estimates of changes in LOS flux near PILs and 
the upward / downward transport of transverse flux. 
Our estimates of flux increase by i7-loop emergence from 
changes in LOS flux near PILs are roughly consistent 
with the net change in LOS flux over the entire sequence. 
Unfortunately, other components of our flux budgeting 
near PILs are inaccurate, in two ways. 

First, our estimates of flux removed by cancellation 
processes (fi-loop submergence and [/-loop emergence) 
are too large (from ~ 10 22 Mx to several times this), 
whether inferred from LOS changes or rates of transport 
of transverse flux. The might be explained if some flux 
is dispersing near PILs that we identify as sites of can- 
cellation, instead of submerging/emerging at these PILs. 
As discussed in Section B~2l identif ying magnetogram fea - 
tures and tracking their evolution (jDeForest et alJ feOO?). 
or imposing additional constraints (e.g., requiring can- 
cellation to persist for multiple time steps, or taking the 
flux lost to be the minimum of that among the two po- 
larities; see IWelsch et al.l 120111 ) could more accurately 
quantify canceled flux. We leave investigation of these 
approaches for future work. 

Second, our estimated rates of transport of transverse 
flux are systematically too large compared to changes in 
LOS flux near PILs, for all processes. As noted in Sec- 
tion 14.21 above, inaccuracies in PIL identification could 
incorporate velocities from siphon flows in our flux trans- 
port rates. In addition, a filling factor / < 1 could also 
play a role, since it would affect |A<I>los| an d |A<I>,jbl| 
differently: while length transverse to the LOS (two fac- 
tors of which enter computation of |A<I>£os|) is scaled 
by y/J ', length along the LOS (which enters computation 
of |A$los|) is not. Hence, if / < 1, one would ex- 
pect |A<&„££| > |A$los|- Given typical discrepancies in 
fluxes by factors of 4 - 6, this would imply filling factors 
of ~ 0.05. This would be implausible in umbrae, where 
the field is expected to be space-filling. We analyze PILs, 
however, which are, almost as a matter of definition, in 
mixed-polarity regions, where the field is expected to be 
more intermittent. The plausibility of a filling factor this 
small along PILs is unclear. 

Another possibility is that flux transport rates are 
overestimated because the transport is not ideal, i.e., 
that the plasma is moving at a different speed than the 
flux transport rate. F or instance, if the flux emergence 
process were not ideal (|Pariat et al.ll2004l) , then the rate 
of flux transport derived from Doppler shifts would un- 
derestimate the true emergence rate, since flux would 
be partially decoupled from the plasma, allowing flux 
to move upward independently of the plasma. While 
magnetic diffusivity would be manifested in that case by 
flux slipping through stationary (or much slower-moving) 
plasma at a rate faster than that suggested by ideal trans- 
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Figure 12. Examples of Doppler shifts along identified PtLs in two HMI LOS magnetograms (separated by about eight hours), both 
without correction for the bias velocity (top row) and corrected using different estimates of vq. Saturation in the color map is set at ±250 
m/s in all panels. Second row: Doppler velocities corrected using vq derived by averaging over all PILs. Third row: Doppler velocities 
resulting by removing the median velocity on all PfL pixels. Fourth row: Doppler velocities corrected using vq derived by averaging over 
high-bias PILs, those with bias fluxes above 3.75 X 10 Mx. Almost all PILs in the top row appear redshiftcd, but red- and blucshifts arc 
more evenly balanced along PILs in the lower rows. 



port, the ideal transport rate is too high in our estimates. 
One possible explanation is that unmagnetized plasma, 
which could be systematically faster-moving than mag- 



netized plasma, could be contributing significantly to 
Doppler signal. 
To quantitatively characterize cases of apparently non- 
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Figure 13. Frequencies of emergence and submergence processes, 
as functions of the number of pixels along all identified PILs at 
all time steps, for three possible corrections to Doppler velocities: 
the average estimates of vq , derived using Btrs with -B t hr = 50 Mx 
cm -2 and d = 5 (top); the median of Doppler velocities on pixels in 
identihed PILs (middle); and the average high-bias estimates of vo, 
again using -Btrs with B t ^ T = 50 Mx cm -2 and d = 5 (bottom). In 
all panels, the black curve shows the distribution of sizes of all PILs. 
Processes involving blucshifts (Q-loop and U-loop emergence) are 
plotted in blue, while processes involving rcdshifts (Q-loop and U - 
loop submergence) are plotted in red. Processes that remove flux 
via cancellation (U-loop emergence and Q-loop submergence) are 
plotted with dashed lines. 

ideal evolution, one could employ a simple non-ideal 
model, assuming a uniform magnetic diffusivity r\ is 
present. This would imply a corresponding non-ideal 
electric field parallel to the PIL (in addition to the ideal 
electric field along the PIL already used to estimate 
A$ vBL ), given by 

c^" = vldwsB^ - d ± B LOS I , (33) 

where, as above, x_l locally points across the PIL (see 
equation [T0|) . For PILs where |A<J> LO s| > \&®vBl\, 



we could attribute the excess A$los to rj 9j_-Blos m ~ 
tegrated along the PIL, taking &losB± to be negli- 
gible, as it would be if there were a current sheet 
along the PIL (e.g., iLitvinenkol 119991 ) . Values of r\ de- 
rived this w ay could be comp a red to the assume d val- 
ues used by iLityinenkd (|1999f ). iChae et all (|2002f l. and 
iLitvinenko et al.l (|2007n . For PILs where |A$ LO s| < 
\A& v bl\, we could attribute the excess A& v bl to over- 
estimation of the transport rate of transverse flux along 
the LOS, due to slippage of transverse flux through the 
moving plasma in the presence of diffusivity. Because 
<9los-B_l is unknown for such cases, however, we could 
not estimate T). 

5.2. Empirical Estimates of Uncertainties 

After filtering azimuths and removing Doppler veloc- 
ities due to solar rotation and convection, we interpo- 
lated the data to sub-pixel accuracy for two purposes: 
(i) image co-registration, to remove frame-to- frame jitter, 
primarily from quantized changes in the cutout window 
boundaries (determined by pixel size) imposed on the 
smooth rotation of the underlying solar features; and (ii) 
Mercator projection of the spherical data onto a uniform 
grid in a Cartesian frame. While interpolation propa- 
gates errors from pathological pixels, and can introduce 
other artifacts (e.g., smoothing) into the data, it should 
nonetheless be undertaken, since significant whole-frame 
shifts were found (some frame-to- frame shifts were on the 
order of 1 pixel, and cumulative shifts showed both sec- 
ular evolution and jumps) in the x and y directions by 
Fourier cross-correlation in the -Blos image sequence. 

Since azimuth filtering and interpolation modify the 
input data, we characterized uncertainties in the dataset 
empirically, instead of relying upon error arrays pro- 
duced in the HMI processing pipeline. To do so, we 
constructed histograms of the LOS and transverse fields, 
-Blos and -B trs respectively, at each time step, and sta- 
tistically characterized these distributions 

In each distribution of -BloSj we interpret the weak- 
field core of distrib ution as a measure of noise (e.g., 
iDeForest et al.ll200~7h . which we quantified by fitting the 
core of the distribution (±6 Mx cm~ 2 ) with a Gaus- 
sian of width ol , and by finding the flux density at half 
max (HM) in the left and right wings. Statistically, the 
mean and standard deviation of fitted ctl's were 4.1 ±0.3 
Mx cm -2 and combined left- and right-widths at HM 
were 4.7 ± 0.5 Mx cm -2 , but fluctuations with a peri- 
odicity around 6 hours were observed. In SDO's geosyn- 
chronous orbit, (i) Doppler motion toward and away from 
the Sun and (ii) thermal cycling due to the Earth's day- 
side albedo would both exhibit 12-hour periodicities, and 
would be a quarter cycle out of phase, so could combine 
to produce variations in instrumental response with such 
a characteristic 6-hour periodicity. This analysis roughly 
agrees with the HMI Team's estimate of the uncertainty 
of 6 Mx cm -2 , and suggests that the uncertainty in the 
LOS field component assumed above, 20 Mx cm -2 , prob- 
ably overestimates the true uncertainty. 

Since B trs is a magnitude, it is necessarily positive, and 
we find the distribution peaked away from zero. If the 

7 We created animations of these histograms 
for the data sequence, which are online at 
http://solarmuri.ssl.berkeley.edu/~welsch/public/data/AR11158/Noise/ 
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distributions of each of the three observed field compo- 
nents were Gaussian, with the same width a, then the 
distribution of transverse field strengths would be that 
of a 3D Maxwell distribution, oc S| ra exp(— B^ rs /2a 2 ). In 
fact, we find that a Gaussian displaced from zero (fit- 
ted to data within ±25 Mx cm -2 of the tallest peak) 
consistently fit each distribution of -Btrs more closely (as 
quantified by \ 2 over the whole distribution) than a 3D 
(or 2D) Maxwellian. It should be noted, however, that 
the shape of the distribution varied substantially over 
the observing sequence — e.g., exhibiting two peaks at 
some time steps — so the Gaussian fits, although better 
than the Maxwellians, were still poor in many cases. The 
means and standard deviations of the Gaussians' fitted 
centroids £?to and widths at were 81 ± 17 Mx cm -2 and 
29±6 Mx cm~ 2 , respectively. Both also exhibited period- 
icities, though with more complicated temporal structure 
than with fits to the distribution of -Blos > with the clear- 
est period at ~ 24 hr. This analysis suggests that our 
assumed uncertainty in the transverse field component, 
90 Mx cm" 2 , is reasonable. 

6. SUMMARY 

We have demonstrated a technique to estimate any 
spatially-constant bias velocity vq present in measured 
photospheric Doppler velocities using a sequence of HMI 
observations of the vector magnetic field near disk cen- 
ter. Such biases can arise from instrumental effects (e.g., 
thermal variations in in strument componen ts) or from 
the convective blueshift (|Dravins et alj|1981l ). Our tech- 
nique uses consistency between changes in LOS flux near 
PILs close to disk center and the transport of transverse 
flux implied by Doppler velocities along those PILs. 

We also noted that an esimate of the bias velocity can 
be made more simply, from the median of Doppler veloci- 
ties along all pixels in LOS PILs close to disk center. This 
approach does not require measurement of the transverse 
field; only an LOS magnetogram and a Dopplergram are 
required. One need only identify PILs, and take the me- 
dian Doppler velocity on all PIL pixels. Since this ap- 
proach ignores our flux-matching constraint in equation 
(fT8"|) . it could be biased by a strong episode of flux emer- 
gence or submergence. 

Accurately calibrated Doppler velocities can be used 
to better unde rstand photos pheric dynamics prior to and 
during CMEs (jSchuckll2010h and to improve estimates of 
the photospheric electric fiel d and thereby esti mate the 
photospheric Poynting flux ()Fisher et al.ll2012T) . and to 
differentiate between possible mechanisms of flux cancel- 
lation (51-loop submergence, U-loop emergence, or recon- 
nective cancellation). 

For the dataset we studied, we found clear evidence of 
bias velocities, though estimates of v ranged from less 
than 100 m s" 1 to ~ 500 m s~ x over the course of the 
run. Only the larger estimates of vq — in the range of 
350 ± 85 m s _1 — were consistent with a majority of 
PILs being sites of emergence, in accordance with the 



substantial increase in total unsigned magnetic flux over 
the duration of our data set. 

Unfortunately, flux budgets based upon our estimates 
of the rates of emergence and submergence are very noisy, 
and likely subject to large errors. Consequently, im- 
proving our approach to constraining the bias velocity, 
or identifying new constraints, would be desirable. We 
plan a careful investigation of systematic sources of un- 
certainty in our method (e.g., better quantitative charac- 
terization of flux cancellation and changes in LOS fluxes 
near PILs), or possible bias (e.g., greater errors in shorter 
vs. longer PILs, or weak-field vs. strong- field PILs). 

It would also be worthwhile to compare results of our 
electromagnetic calibration with calibrations based upon 
fitting the center-to-limb component of full-disk HMI 
measurements of Doppler veloc ities (iSnoderassI 119841 : 
lHathawavi[l99l lHathawav etalll255ot lSchuckH2010[ TAs 
noted above, center-to-limb fitting has a major physical 
uncertainty: the physics of the center-to-limb variation 
in Doppler shift in the particular spectral line used by 
HMI (or any other spectral line) involves detailed inter- 
actions between height of formation, the height of con- 
vective turnover, the variation with viewing angle of the 
average convective flow speed , and the variation with 
viewing angle of optical depth ( )Carlsson et al.ll2004[ ) . Di- 
verging flows tangent to the photosphere in granules can, 
depending upon optical depth in granules at the forma- 
tion height of the line, produce a convective blueshift 
at the limb, because diverging flows on the near sides of 
granules approach the viewer, while receding flows on the 
far sides of the granules are at least partially obscured by 
the optical depth of the granules (see Figure [5]). Hence, 
finding the center-to-limb shift in a line does not imply 
that the line does not still possess an overall additional 
shift. 

We note that, in addition to the Dopplergram data an- 
alyzed here, the HMI Team also produces independent 
estimates of Doppler velocities in the Milne-Eddington 
inversion process. While comparisons between ME and 
Dopplergram velocities show them to be quite similar 
(Yang Liu, private communication 2011), these indepen- 
dent measurements can provide additional information, 
e.g., better characterization of uncertainties. 

Data and images are courtesy of NASA/SDO and the 
HMI science teams. We thank Stanford University's 
Yang Liu, Jesper Schou, and Sebastien Couvidat, and 
NWRA's K.D. Leka for readily answering queries about 
HMI data processing. SDO/HMI is a joint effort of many 
teams and individuals, whose efforts are greatly appre- 
ciated. This research was funded by the NASA He- 
liophysics Theory Program (grants NNX08AI56G and 
NNX11AJ65G), the NASA Living- With-a-Star TR&T 
Program (grant NNX08AQ30G), and support from NSFs 
National Space Weather Program (grant AGS-1024862). 
The authors are grateful to US taxpayers for providing 
the funds necessary to perform this work. 



APPENDIX 

GEOMETRIC EFFECTS AWAY FROM DISK CENTER 

We have asserted that the calibration method we describe here is only valid near disk center. To better understand 
why, we consider a particular configuration of emerging magnetic flux in Figure [T4J where it is illustrated in two 
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different orientations with respect to the viewer. This configuration is highly idealized: we assume the configuration 
is invariant in the direction into the page, and that emerging held lines have the same shape as held lines that have 
already emerged. In both panels, the viewer is assumed to look down from above the top of the figure. 

In the upper left panel, a plane tangent to the photosphere at the radial- held PIL is normal to the LOS. At this 
particular viewing angle, the photospheric normal direction and the LOS coincide. Away from the PIL, the held is 
uniformly tilted from the LOS by the inclination angle O = 45°. In the upper right panel, the normal to the plane 
tangent to the photosphere at the radial-field PIL is tilted with respect to the LOS by an angle 9 — 45°, which is also 
the tilt of the legs of the emerging flux system with respect to the photosphere. 

We assume the flux emergence is ideal, such that the electric held obeys equation (|3|) , and flux is therefore frozen to 
the plasma. The rate of flux emergence by a radial flow can therefore be inferred from observations of the upper-left 
panel in two ways: (1) from the increase in the the total unsigned LOS flux (the integrated absolute value of the curve 
plotted in the lower left panel); and (2) from the product of the transverse field at the PIL and the Doppler signal at 
the PIL. The equivalence of these measures of the rate ideal flux emergence underlies the Doppler calibration method 
presented in this paper: we assume any inconsistency between (1) and (2) is due to a bias in the Doppler zero point. 

In contrast, flux emergence can only be unambiguously inferred by observations of the upper right panel one way: 
from the increase in the the total unsigned flux (either LOS flux or radial flux). This is because the magnetic field at 
the radial-field PIL has a component along the LOS, implying Doppler shifts at this PIL might contain a contribution 
from flows parallel to B there. Flows parallel to B, however, do not drive flux emergence, and cannot be distinguished 
from Doppler shifts due to flows perpendicular to B. 

It is also clear that the relative orientation of the negative-polarity flux and the LOS implies that an observer would 
not see an increase in negative-polarity LOS flux as emergence proceeds. Hence, observations of the emergence of 
tilted flux in LOS magnctograms can exhibit strong flux imbalance. Measurements of changes in the radial flux, from 
vector magnetograms, would be balanced in this idealized case. 
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Figure 14. Upper left and right panels depict a highly idealized emerging flux configuration, assumed to be invariant in the direction into 
the page, viewed in two different orientations with respect to an observer, who is assumed to be looking down from above the top of this 
figure. In the upper left panel, measurements of Doppler velocities on the radial-held PIL and transverse Gelds at the PIL can be used to 
infer the rate of flux emergence. In the upper right panel, the magnetic Geld at the radial-Geld PIL has a component along the LOS, so Gows 
parallel to B (which do not drive Gux emergence) can contribute to Doppler shifts, implying the Doppler velocity cannot unambiguously 
determine the rate of Gux emergence at this PIL. In both cases, Gux emergence can be detected from increases in total unsigned LOS Gux 
(the curves plotted in the bottom left and right panels), and quantitatively estimated from changes in the total unsigned radial magnetic 
Gux. 
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